Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



O 

o 



o 



> 
o 

t> 
in 
o 
\o 
o ■ 

Or 

O ■ 

a: 



13 
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ABSTRACT 

New measurements of the z ~ 6 intergalactic C iv abundance are presented, using moderate reso- 
lution IR spectra of two QSOs taken with GNIRS on Gemini South. These data were systematically 
searched for high redshift C iv absorption lines, using objective selection criteria. Comprehensive 
tests were performed to quantify sample incompleteness, as well as the rate of false positive C IV iden- 
tifications. The trend of constant fic iv( z ) observed at z ~ 2 — 5 appears to continue to z ~ 6, the 
highest observed redshift. The C iv sample is also con s istent with the redshift-invariant form of the 
C IV column density distribution reported bv iSongaTIal l)2001j) at lower redshift, although with fairly 
large uncertainties due to a smaller sample size and noisier infrared data. The constant value of IV 
does not necessarily imply that the IGM was infused with an early metallicity "floor," but the pres- 
ence of early C IV does indicate that heavy-element enrichment began < 1 Gyr after the Big Bang. 
The lack of a decline in f2c iv at high redshift may indicate that integrated C iv measurements are 
sensitive to the instantaneous rate of feedback from galaxy formation at each epoch. Alternatively, 
it could result from a balance in the evolution of the intergalactic gas density, ionization conditions, 
and heavy-element abundance over time. 

Subject headings: 



1. INTRODUCTION 

Recent CMB measurements of the electron scattering 
optical depth suggest that star formation first began to 
ionize the IGM at z » 11, a bout 425 Myr after the 
Big Bang j|Spergel et al.ll2006j) . The subsequent epoch 
be tween z = 11 a nd the end of reionization at z 
6 llFaxi et al. 1 120021 
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120011) encompasses only ~ 500 Myr. This interval be- 
gins to approach the dynamical timescale for massive 
galaxy formation and supernova feedback, so the galax- 
ies formed at z ~ 6 — 10 may be among the first objects 
to produce the heavy elements seen in the IGM at z ~ 3 
ijSchave et alJl2053 ISimcoe et ai]l200|) . 

If the large scale enrichment of the IGM actually be- 
gan within this short window, one would expect to see 
either a significant downward trend in heavy-element 
abundances, or at least an increase in abundance scatter 
approaching z ~ 5 — 6. Even if the heavy elements are 
removed from early galaxies by strong winds, there are 
constraints on the distance that these winds will cover 
based simply on travel time (e.g. D ~ 10 kpc per lOOMyr 
for a 100 km s -1 wind). It is also likely that early galactic 
winds would stall at R < 100 kpc because of mass load- 
ing from both the ISM and the IGM, which is 300-1000 
times more dense at these redshifts than in the present 
day. 

1 Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership: the National Science Foun- 
dation (United States), the Particle Physics and Astronomy Re- 
search Gouncil (United Kingdom), the National Research Coun- 
cil (Canada), CONIC YT (Chile), the Australian Research Council 
(Australia), CNPq (Brazil) and CONICET (Argentina) 

2 MIT Center for Space Research, 77 Massachusetts Ave. #37- 
673, Cambridge, MA 02139, USA; simcoe@mit.edu 

3 Pappalardo Fellow in Physics 



To date, no abundance decline has been observed in 
the IGM towards high redshifts. Over the z ~ 2 — 4 
range where one can estima te ionization-correct ed car- 
bon abundances from C IV , iSchave et al.l l)2003|) found 
no significant redshift evolution i n [C/H] across a sam- 
ple of several sightlines. Likewise, iSongailal l)20Qll l2005[) 
and lPettini et alJ l)2003|) found no evolution in iv, the 
integrated contribution of C iv gas to closure density, 
between z = 1.8 and z = 5. The latter result is partic- 
ularly surprising since flc rv depends on both the over- 
all carbon abundance and the C iv ionization fraction 
nc iv/n-c- Naively one would expect both of these quan- 
tities to vary over such a broad redshift range, because 
of ongoing chemical enrichment and the evolution of the 
ionizing background radiation spectrum. 

Despite the discov ery of severa l z > 6 QSO s in the last 
few years l|Fan et al.ll200ll 120031 l200l. |2006|) . intergalac- 
tic absorption-line studies have been limited to z < 5 
because at higher redshift the C iv doublet moves into 
the near-infrared. Absorption line spectroscopy is much 
more challenging in this regime because of increased de- 
tector noise, OH emission from the sky, and more severe 
telluric absorption. These foreground problems can be 
mitigated by employing higher spectral resolution, but 
the associated penalties in sensitivity require long inte- 
grations to overcome. 

In this paper I present J band GNIRS spectra of 
two z ~ 6 QSOs at moderately high spectral resolu- 
tion (Av = 60 km/s). These spectra are searched for 
C iv absorption systems at 5.4 < z < 6.2, and the results 
used to extend existing measurements of fic iv to z ~ 6. 
In section |5J I describe the observations and reduction 
methods. Section |2| details the C iv search procedure 
and simulations to quantify the effects of incompleteness 
and/or spurious C IV detections in the data. Sections 0] 
and present the results for f2c iv and a discussion of 
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Fig. 1— C iv spectra of z ~ 6 QSOs SDSS1306+0356 
(z em = 6.002) and SDSS1030+0524 (z em = 6.272). Solid line 
shows continuum estimate, and light green line shows la error es- 
timates, including the effects of Poisson noise from sky emission 
and telluric absorption. Night Sky emission is shown in bottom 
panels, and the C IV search ranges are indicated with labeled hor- 
izontal bars. The redshift of the C IV emission peak is shown with 
a vertical dash. 



their implications and limitations. 

2. OBSERVATIONS 

Figure n displays the C iv spectral regions for the high 
redshift QSOs SDSS1306+0356 (z cm = 6.002 as mea- 
sured from the C IV emission peak) and SDSS1030+0524 
( z cm = 6.27 2). The data were obtained with the GNIRS 
instrument ijElias et al.lll998j) . operated in queue mode 
on Gemini-South during semesters 2005A and 2006A. 
The instrument was configured with the cross-dispersion 
prism, the 110 line/mm grating, and a 0.45 arcsecond 
slit. Two grating positions were observed for each QSO 
to provide complete spectral coverage up to the emission 
wavelength of C iv , although the signal-to-noise ratio 
(SNR) is not completely uniform between setups due to 
varying conditions on different nights within the queue. 
Each QSO was observed for approximately 8 hours per 
setup. 

The observations were split into sequences of 4 x 10 



minute integrations, nodded four arcseconds along the 
slit in an ABBA configuration. The data were reduced 
using customized IDL software developed by the author. 
Individual frames were filtered for detector pattern noise 
and flatfielded, and the ABBA frames were combined to 
remove signatures of the sky, dark current, and charge 
persistence. Each 2D pixel was assigned a wavelength by 
fitting a second order polynomial to the centroids of OH 
sky lines , which were compared with the line lists com- 
piled bv iRousselot et ahl 11200 0). Residual errors in the 
sky subtraction were corrected by a second-pass fit of a 
b-spline function to the subsampled grid of wavelength 
versus pixel intensity. Objects were identified and ex- 
tracted using an optimal Gaussian profile, with errors de- 
termined by direct calculation of the sky RMS outside of 
the object profile (the errors were not dominated by shot 
noise from the object). The data were flux calibrated 
and corrected for telluric absorption using observations 
of standard stars obtained throughout the night at sim- 
ilar airmass as the QSOs. Finally, the extracted spectra 
were scaled to match in flux, and added together on a 
ID grid of pixels with constant velocity spacing of 25.3 
km/s. Exposures from different nights were weighted ac- 
cording to their (SNR) 2 , and the wavelengths converted 
to vacuum heliocentric coordinates. The final resolution 
of the combined spectra was measured as R = 5000 or 
60 km/s FWHM (2.35 pixels) using Gaussian fits of sky 
lines in both the reduced ID and the raw 2D frames. 

3. ANALYSIS 

The data in Figure ^ are among the first z > 5 IR 
QSO spectra with resolution suitable for absorption line 
measurements; however, they are still a factor of ~ 10 
coarser than high resolution optical spectra at z ~ 3 
and have lower SNR by a factor of several. At Av = 60 
km/s per resolution element, they are sufficient to resolve 
the C iv doublet (with characteristic Av ~ 500 km/s). 
However, individual C iv components at lower redshift 
have typical intrins ic velocity dispersions of b ~ 10 km/s 
ijRauch et al.l lT996). and hence will generally not be re- 
solved in these spectra. Since the line broadening is dom- 
inated by instrumental effects, I have characterized the 
absorption systems through a curve-of-growth analysis 
with attention to possible saturation effects, rather than 
by direct Voigt profile fitting. However, all of the results 
were compared against Voigt profiles generated using in- 
dependent software as a final consistency check. 

The data were divided into two separate redshift inter- 
vals for the analysis, with the interval boundaries deter- 
mined by the SNR of the data. The redshift ranges are 
indicated in Figure^by horizontal bars. In both objects 
an increase in the pixel noise level is seen at the transi- 
tion between two spectral orders at A ~ 1.04 /im. This 
established the lower boundary of the high-redshift bin 
at z = 5.72. The upper boundary of the high-redshift 
bin was set at z = 6.21, or 2500 km s _1 below the emis- 
sion redshift of SDSS1030+0524. It is common prac- 
tice in absorption line surveys to omit this last portion 
of spectral coverage because any absorbers there would 
be affected by local radiation from the QSO. For this 
same reason the sarch range of SDSS1306+0356 was lim- 
ited to z = 5.72 — 5.93 for the high redshift bin. Nei- 
ther SDSS1030+0524 or SDSS1306+0356 contains any 
C IV lines within 2500 km s _1 of its emission redshift. so 
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the results are not sensitive to this choice. Inclusion of 
the extra spectrum would amount toa~ 15% (0.06 dex) 
increase in pathlength. 

Because of its lower emission redshift, SDSS1306 does 
not cover the entire range of the high redshift bin. Hence 
the effective (pathlength weighted) central redshift of the 
bin falls at z = 5.92, slightly below the nominal midpoint 
of z = 5.97. The low redshift bin ranges from z — 5.40 — 
5.72, with its lower bound chosen to avoid contamination 
from interloping Mg II at A « 99001 in SDSS1306+0356. 
Both objects cover its full width. 

Because of the IR data's lower spectral resolution, 
lower SNR, and shorter pathlength relative to z ~ 3 opti- 
cal surveys, this initial abundance measurement at z ~ 6 
will not be sensitive to weak C IV systems in the very ten- 
uous IGM. However, the data quality is high enough to 
detect the stronger C iv systems commonly seen in lower 
redshift data with high confidence. Because the slope 
of the C IV column density distribution is flatter than 
iv (Songaila 2001), it is these sys tems which provide 
the dominant cont ribution to flc iv ijPettini et al.ll2003l 
iSimcoe et alJ 120061. and this is the main measurement 
presented in Section l4~2l 

3.1. Identification and Measurement of C IV Lines 

Each of the QSO spectra was searched systematically 
for C IV doublets using an automated software routine. 
First, an estimate of the continuum was constructed by 
interpolating an interactively fit cubic spline. The mini- 
mum spacing of the spline knots was set at 500 km s _1 to 
avoid overfitting the continuum on smaller scales. 

The entire redshift range was then scanned for ab- 
sorption feaures with > 3c significance in equivalent 
width (EW) over a resolution clement. Each flagged 
feature was considered a possible 1548. 2A component of 
the C iv doublet, and the EW was then calculated for 
the corresponding 1550. 8A region. If the 1550. 8A EW 
was also > 2<j, the system was flagged as a candidate 
C iv doublet. The list of candidate doublets was then 
examined by hand to discard systems that were obvi- 
ously contaminated by sky subtraction residuals or noise, 
or which had doublet ratios inconsistent with C iv (i.e. 

Wi55i > Wl548)- 

Application of this procedure resulted in nine C IV de- 
tections between the two sightlines, ranging from z = 
5.451 to z = 6.175. For each of these systems, I esti- 
mated a more accurate rest EW of C iv (Wr,c iv) by 
fitting a Gaussian absorption profile to the data, with 
its width fixed at one resolution element. Plots of each 
system and its best-fit absorption profile are shown in 
Figures through 0] with absorption line redshifts and 
EWs listed in Tables 1 and 2. 

Some of the high redshift C iv lines are highly signif- 
icant detections, while some are more marginal. Along 
the SDSS1306 sightline (which has a higher SNR) only 
three systems are present. The system at z — 5.674 
is a strong detection, although it is blended with noise 
features in the data. The systems at z = 5.451 and 
z = 5.927 are much weaker and even with a high SNR 
over these portions of the spectrum, the C iv 1550. 8A de- 
tection significance is low. 

Toward SDSS1030+0524 the search algorithm identi- 
fied six C iv systems. Two strong detections are found 
at z = 5.738 ± 0.005 and z = 5.827 ± 0.005. The sys- 
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Fig. 2. — C IV detections along the line of sight towards 
SDSS13064-0356. Best-fit absorption models are shown with solid 
lines. The upper two systems are contained in the low-redshift 
sample; the bottom system is a marginal detection, and it is the 
only C IV line identified at high redshift toward this QSO. 



tern at z = 5.738 contains two significant components, 
with the 1548. 2A component of one line blended with the 
1550. 8A component of the other. It is found in a noisy re- 
gion of the spectrum, where systematic errors in the sky 
subtraction are non-negligible. The z = 5.827 system 
contains three highly significant doublets. This complex 
also contains another strong line to the red which could 
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Fig. 3. — C IV detections along the line of sight towards 
SDSS1030+0526. Top system at z = 5.458 is in the low redshift 
sample, all others are in the high redshift sample. Best-fit absorp- 
tion models are shown with solid lines. 
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Fig. 4. — C IV detections along the line of sight towards 
SDSS1030+0526 (Continued from previous figure). 



not be identified with any ion commonly seen in QSO 
spectra at any redshift. All well-known QSO lines with 
rest wavelengths between 1450A (corresponding to the 
QSO's rest frame) and 6000A (i.e. past the NaD doublet) 
were considered without finding a match. The presence 
of this absorption was verified in the 2D trace over multi- 
ple nights, and checked to rule out artifacts from telluric 
absorption or features in the flat field. It was also checked 
against a high-resolution optical spectrum of SDSS1030 
(kindly provided by G. Becker) without resulting in any 
cross-identification. The most likely interpretation ap- 
pears to be a mixture of three C iv components with one 
additional unidentified line. 

Only two weak C iv systems are identified at z > 6. 
They are each near the detection limit of the data, and 
both are found in relatively clean and sky-free portions of 
the spectra. It is intriguing that the strong C iv systems 
generally seem to be found at z < 5.9, while at z > 5.9 
(where the SNR of our data is highest) the C iv sys- 
tems are weaker. While this could be consistent with 
some evolution of the population, this effect could also 
be caused by shot noise along a short redshift path, as 
described in more detail in Section l4~T1 

Some, but not all, of the C iv doublets show evidence of 
instrumentally unresolved saturation in the physical ab- 
sorption profile. Figure shows the C iv curve-of-growth 
(COG) of rest EW (W r ) versus Nq iv ; the bottom panel 
displays the EW ratio between the two components as a 
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Fig. 5. — Curve-of-growth calculations for C IV . Top panel shows 
the relation of W r to Nq iv for both components of the doublet. 
Bottom panels shows the equivalent width ratio of the two doublet 
components as a function of Nq iy, for three typcal choices of b 
parameter. 



function of Ac iv for various values of the C iv b param- 
eter. For reasonable values of bciv, the 1548. 2A com- 
ponent begins to saturate near Ac iv ~ 10 13 cm~ 2 , and 
the doublet ratio approaches unity as the 1550. 8A com- 
ponent also begins to saturate at Ac iv ~ 10 14 cm -2 . 

Column 7 in Tables ^ and [21 shows the measured dou- 
blet ratio for each C iv detection. Because of the data's 
sensitivity limitations, most of the detected systems lie 
in the non-linear regime between 10 13 < Ac iv ^ 10 14 
cm~ 2 , where the Ac iv — W r mapping is degenerate with 
the C IV b parameter. However, even for some of the 
stronger systems in the data with W r ~ 200 — 300 
mA (e.g. the complex at z — 5.823 — 5.830), the dou- 
blet ratios provide useful upper bounds on Nc iv Often 
the ratio is larger than one would guess from measur- 
ing W r (1548) alone, and in these cases it is likely that 
the observed C iv line is actually an unresolved blend of 
several narrow components with more complex velocity 
structure. 

At z ~ 3, the typical C iv line has a velocity dispersion 
of b 10 km s _1 , and over 80% of C IV components 
have &c iv < 16 km s -1 l|Rauch et al.l[l996j) . Assuming 
be iv = 10 km s _1 for our systems, the doublet ratio 
can be used to estimate Nciv- Then, the number of 
unresolved components can be inferred by comparing the 
measured value of W r with the single-component value 
from the COG. The rightmost columns in Tablesn an d|21 
show the Nc iv value deduced from the EW ratios, along 
with the corresponding number of components required 
to match the total EW for each system. 

The Nq iv values for individual components range be- 
tween 10 1315 and 10 143 . Different choices of 6c iv can 
change specific values of Nc iv by up to ~ 0.2 dex; ran- 



dom measurement errors in the EW could add another 
0.2 dex of uncertainty, though these errors would be un- 
correlated. However, the total column density of each 
C IV system is fairly insensitive to the specific choice of 
b. A larger b parameter leads to fewer unresolved com- 
ponents with higher Nc iy, while a small b leads to more 
components with smaller column densities, but the sum 
of the columns matches to within < 0.1 dex for any choice 
of b between 5 and 20 km s _1 . 

3.2. Estimates of Sample Completeness 

IR absorption line searches require particular attention 
to characterize sample incompleteness and false-positive 
identifications. In addition to normal SNR variations in 
the data, one encounters many spectral regions where 
Poisson noise from OH sky emission or residuals from 
imperfect telluric absorption correction render pixels un- 
usable for cosmological absorption measurements. At 
R = 5000, the GNIRS spectra resolve sky features well 
enough to work in many of the dark regions between fe- 
tures from the night sky. 

To quantify the completeness of the z ~ 6 C iv sam- 
ple as a function of Nc rv, we have performed a series of 
Monte Carlo simulations. For each redshift bin in each 
sightline, a sample of artifical C iv doublets was created 
with Nc iv ranging from 10 12 — 10 16 cm~ 2 and redshifts 
selected at random. Each artifical line was added individ- 
ually into the real spectrum, and the data were searched 
using the same EW criteria to see if the artificial line was 
recovered. One thousand trials were performed for each 
column density. 

Figure displays / comp , the fraction of lines recov- 
ered as a function of Ac iv for each of the sightlines. 
For most of the redshift path, the data are largely com- 
plete (> 70%) at A c iv > 10 14 cm" 2 . The complete- 
ness drops off rapidly towards Ac i V < 10 13 cm -2 , at 
which point only 5 — 10% of systems are detected. In the 
IR, unlike the optical, it is nearly impossible to achieve 
100% completeness even at thie highest column densi- 
ties because of chance coincidences with OH sky lines or 
strong telluric features. However, for the C IV lines at 
Ac iv ^ 10 14 cm~ 2 , it appears that at R ~ 5000 one re- 
solves night sky foregrounds out of ~ 80% of the total 
available pathlength. 

Sample incompleteness effectively reduces the absorp- 
tion pathlength over which lines of a given column den- 
sity can be detected in the survey. In the discussion 
below, we explicity state when completeness corrections 
are applied to the data. These are implemented by scal- 
ing the relevant pathlengths AA by the factors shown in 
Figure 

3.3. Incidence of False Positive Identifications 

Several of the C iv detections in Tables ^ and |3 ap- 
proach the sensitivity limit of the data, and at these 
column densities one must also consider the possibility 
of false positive identifications from chance noise fluctu- 
ations. This process is much more difficult to simulate 
since in true spectra the noise properties are often not 
perfectly approximated by a Poisson distribution, espe- 
cially in the very instances where correlated pixels lead 
to false detections. 

One method for quantifying the false positive rate is 
to scan the spectra for absorption "doublets" identical to 
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Fig. 6. — Estimates of sample completeness as a function of 
Nc iv i determined from Monte Carlo simulations performed on 
the data. Top panel shows completeness for SDSS1306 over both 
sample redshift bins; bottom panel shows the same for SDSS1030. 
C IV lines are recovered more completely in the higher redshift bins 
because of improved SNR in these spectral regions. 



C iv in every sense, except that the oscillator strengths 
of the two transitions are reversed. Systems identified 
this way will be subject to the same noise properties and 
line spacing as true C IV, but can only be caused by 
correlated noise fluctuations in the data since a doublet 
with W r (1551) > W r (1548) is unphysical. 

I searched for this type of system in the two QSO spec- 
tra using the same EW criteria as for the true C iv search. 
For SDSS1306, the high redshift sample at z = 5.72-5.93 
yielded one weak "false" system. For SDSS1030, there 
was one (rather strong) flase system at z = 5.722, and 
two other weak systems at z = 6.005 and 6.155. Over the 
full z = 5.72 — 6.21 redshift path, the total column den- 
sity of all "true" C iv systems is 10 14 8 cm" 2 , whereas 
the the total for the "false" systems was 10 143 cm -2 . 
Contamination from false positives may therefore com- 
prise up to < 30% of the reported C IV in this bin, but 
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Fig. 7. — The C IV column density distribution function, de- 
fined as f(Nciv) = dM / dNc ivdX , over three different redshift 
intervals. Points are shown with incompleteness corrections; the 
uncorrected values are shown with open circles. Left and center 
panels show the measured CDDF for the nominally defined red- 
shift bins. The right panel illustrates limits for a possible high 
redshift bin, where more absorption pathlength would be required 
to obtain an accurate measurement. 



spurious detections do not appear to dominate over the 
true C iv signal. In the z = 5.40 — 5.72 bin the false 
positive signal is likewise estimated at < 25%. 

If one restricts the search to the z > 6 region, the con- 
tribution from contaminants is no longer subdominant. 
The integrated C iv column for the three "true" C IV sys- 
tems is ^2 Nq iv = lO 13 6 cm" 2 , whereas the two "false" 
systems have Nc iv = 10 137 . This result is consistent 
with the interpretation that all of our z > 6 C iv detec- 
tions are spurious, or alternatively that there is no statis- 
tically significant signal from C iv over z = 6.0 — 6.21 in 
SDSS1030. This is an intriguing possibility, but the ab- 
sorption pathlength is fairly short for this one sightlinc, 
so it could just be the result of statistical fluctuations in 
a small sample. I return to this point in the discussion; 
for the purpose of calculations below I continue to treat 
the z > 6 C iv detections as real, keeping in mind that 
because of contamination they will lead to upper limits 
on the values of any derived quantities. 

4. RESULTS 

4.1. z ~ 6 C IV Column Density Distribution 

A key resul t from lower redshift C iv surveys (e.g. 
ISongailal 12001) is that the column density distribution 
function (CDDF) of C iv — defined as the number of ab- 
sorbers per unit Nq iv per unit absorption pathlength — is 
nearly invariant between z ~ 2 and z ~ 5. Even with 
the small sample available at z ~ 6, one can construct 
a coarse CDDF for comparison with this lower redshift 
data. 

For calculating the absorption pathlength, I use the 
convention that AX = X[z-x) — X(zi), where X(z) is 
defined as 

X ( z ) = iJ—VVm(1 + z)3 + {I a (1) 

with fijvf = 0.3 and Q\ = 0.7. Results from the literature 
are often quoted using the AX for an Einstein-deSitter 
universe; where appropriate I have rescaled the literature 
values by to conform with more recent concordance 

models. 
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Given the small size of the z ~ 6 line sample, I only es- 
timate the CDDF in two Nc iv bins — one ranging from 
log(iVc iv) = 13.00 - 13.75 and one from log(iV c IV ) = 
13.75 — 1 4.5. For consis tency, I have followed the ap- 
proach of HSflngajlaJ |2fJQl(), summing the Nc iv values of 
individual components in each C iv complex to obtain a 
single, total column density for the system. For example, 
the lines spanning z = 5.82374— 5.83012 were treated as 
a single system with New = 10 14 ' 42 , corresponding to 
the sum of the individual components listed in Tabled 

The resulting distributions are shown in Figure and 
Table □ The solid line shows the fit to the C IV CDDF 
derived from the z ~ 3 sample of lSongailal (|2001fl . after 
rescaling to SIm — 0.3. The open circles indicate the 
2^6 CDDF values obtained before correcting for sam- 
ple incompleteness (as described in Section l3"T2"|) . Com- 
pleteness corrected points are shown with solid dots; the 
corrections were made by rescaling AX by the recovery 
fractions shown in Figure El The horizontal error bars 
indicate bin sizes, and the vertical bars indicate Poisson 
counting errors. 

Despite the modest size of the z ~ 6 sample, the data 
do appear to be consistent with a non-evolving form of 
the CDDF at z — 5 — 6. In the primary sample with 
highest data quality at z = 5.72 — 6.21, the data fall di- 
rectly on top of the low redshift fit. In the lower redshift 
bin the SNR is also lower so the sample is substantially 
incomplete (e.g. there is only one system in the the low 
column density bin, and we expect sample contamina- 
tion). However, within substantial errors even this sam- 
ple is consistent with a non-evolving distribution. 

It would be stretching the data somewhat to consider 
these calculations a detailed measurement of the CDDF 
at z ~ 6; a full and accurate measurement would require 
both more sightlines to build up AX at Nc rv ^ 10 14 
cm~ 2 (where our sample is mostly complete), and higher 
SNR data to improve the sample completeness for weaker 
lines at Nqiv ^ 10 13 5 cm~ 2 . However, the fact that we 
detect any lines at all argues against a strong downward 
trend in the C IV abundance between redshift 5 and 6, 
and the similar overall normalization (about one system 
per unit absorption path, per logarithmic unit column 
density) suggests that these is no statistically significant 
evidence for evolution. 

It is notable that within the z = 5.72 — 6.21 bin, most 
of the C IV absorption is concentrated at the bluest end; 
in fact in Section EPl it was suggested that for z > 5.9 the 
data may be statistically consistent with no C IV signal. 
It is tempting to interpret this as a hint of decline in the 
C iv abundance, but there are dangers in this sort of a 
posteriori definition of a redshift boundary, particularly 
for small sample sizes (recall that the redshift bins for 
the main sample were bounded only by random factors: 
a change in the noise properties of the data, and the QSO 
emission redshifts). 

Nevertheless, for argument's sake one can ask whether 
the lack of C iv at z > 5.9 should be surprising 
given the amount of absorption path covered in the 
SDSS1030+0526 spectrum. In the rightmost panel of 
Figure we separate out the z = 5.9 — 6.28 red- 
shift range, and indicate the upper limits derived from 
assuming that there is < 1 system at log(iVciv) = 
13.75 — 14.50 (we detected none) and < 4 systems at 
\og(Nc iv) = 13.0 — 13.75 (we detected 4 systems, but 
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Fig. 8. — Evolution of the C IV contribution to closure density. 
New points are shown with round dots at z = 5.40 — 5.72 and z = 
5.72 — 6.21. Verical error bars r epresent 90% con fidence intervals. 
Results at z = 1.8 — 5.5 from Songailal I2001T are shown with 
squares, and the z = 4.5 — 5.0 point from Pctti iri et alJ J2003) is 
shown with an open triangle. All values are scaled to h = 0.71, 
Q M = 0.3, and Q A = 0.7. 



also 3 false positives). Taking just the strong sys- 
tems at N c iv = 10 13 - 75 ~ 14 - 5 which should be easily 
detected, the completeness-corrected absorption path- 
length in the SDSS1030 sightline is only marginally large 
enough (AX — 1.16) to expect a single strong C iv de- 
tection, so the absence of strong C iv could be the result 
of a statistical fluctuation. The probability of observ- 
ing zero strong systems in SDSS1030 at z > 5.9 can be 
estimated from f(N C iv) (the C IV CDDF) as: 

P(0) = exp [-/(JV C iv)AXAJV c IV ] « 50%. (2) 

So, for the present, the possibility of a decline in 
C iv abundances at z > 5.9 remains intruiging but the 
sample is too small to make conclusive statements. If the 
CDDF remains constant at z > 6, then there is a ~ 75% 
chance that a strong system will be detected in the next 
sightline to be searched at z > 6. 

4.2. riciv' Contribution to Closure Density 

The integral of the CDDF ( weighted by Nc rv ) can 
be used to estimate the volume-averaged density of 
C iv atoms over a specified redshift interval. This quan- 
tity is customarily normalized by the critical density p c 
to yield f2c rv; which is calculated from observed quanti- 
ties as: 

1 E^Civ 

"Civ = — me iv c ^ T v 

pc iEAi 



(3) 



where mc iv is the mass of the C iv ion. This integrated 
C iv density provides a conveient metric for study i ng the 
redshift evolution of C iv abundances. iSongailal 1)200 If) 
has com piled measurements of fic rv i n the z = 2 — 5 
interval; iPettini et alJ l)2003t) have also published a mea- 
surement at z = 4.5 — 5.0. Figure |H1 displays these prior 
results, for Hq = 71 km s^Mpc -1 , and AX rescaled as 
above for VLm = 0.3 and f2^ = 0.7. 

The rightmost points show the two new estimates of 
O c IV at z = 5.40 - 5.72 and z = 5.72 = 6.21 from this 
analysis. Horizontal error bars indicate the bin ranges, 



and vertical error b ars represent 90% confidence inter- 
vals calculated as in IPettini et alJlfc OO.U using the error 
estimator of lStorrie-Lombardi et alJ l)1996|) : 



C IV J 



n 



(90%) = 1.64 x 



C iv 
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The new points in Figure |S] have been corrected for sam- 
ple incompleteness by scaling each line's column density 
in the sum from Equation[3]by the pathlcngth over which 
it could be detected. It is also corrected for spurious con- 
tamination by subtracting off the total column density of 
false-positive lines identified in each bin (see Section lST^|) . 
These two corrections are both of order ~ 0.2 dex but 
act in opposite directions. 

The z ~ 6 points follow the trend of a constant 
ficiv established at lower redshift, and do not show 
evidence for a decline at z > 5. While some hints of 
a declining Qc iv were present in the optical data at 
z ~ 4.5 — 5.5, this may be an artifact of shot noise, 
as was pointed out by Songaila in the work where it 
was reported. The z = 5 — 5.5 range is still particu- 
larly difficult for such measurements because of the low 
sensitivities of both CCDs and IR focal plane arrays 
at A ~ 0.98/im, and the substantial atmospheric fore- 
grounds. Both Songaila's measurement at z = 5 — 5.5 
and the new point at z = 5.4 — 5.72 are limited by shot 
noise. A pathlength-weighted average of the two points 
yields fi c i V (5.00 - 5.72) = 3.6±|;| x 10~ 8 , consitent with 
the full redshift range, though still with substantial error 
bars. In the present sample the z — 5.72 — 6.21 point 
is more reliable, since it incorporates more absorption 
path, it is situated in a cleaner portion of the Y band, 
and it represents data with a higher SNR. The error bars 
for both z ~ 6 points are larger than for the lower red- 
shift data because we have only sampled two sightlines; 
however since redshift path is a comoving quantity pro- 
portional to (1 + z)dz the effective path per sightlinc 
increases at higher reshift for Qc iv measurements. 

Given the resolution and SNR of the z ~ 6 spectra, 
even a single C iv detection is enough to bring Qq i V up 
to the ~ 10 -8 level. If the C iv density was in rapid 
decline at z > 5, no C iv lines should have been detected 
in the GNIRS spectra, and in this case our sensitivity 
would have been sufficient to set upper limits roughly 
a factor of 3 below the ~ 3 x 10~ 8 level seen at lower 
redshift. Since several C iv lines are detected, even if a 
fraction of these are false-positives, it seems unlikely that 
they could mask an order-of-magnitude decline in Qq iv 
at z ~ 6. 

5. DISCUSSION 
5.1. Implications for Cosmic Abundances 

The z ~ 6 C iv systems detected in SDSS1030+0524 
and SDSS1306+0256 yield an fl c rv value of « 3 x 10~ 8 , 
which is essentially identical to the value observed at 
lower redshift. However, this measurement reflects only 
the C IV abundance, and does not take into account the 
fraction of carbon atoms in other ionization states. 

A correction for ionization requires estimates of fc rv = 

iv 

/n c — a factor which varies with redshift and envi- 
ronment and can be large. In ionization equilibrium, the 
exact value of fc iv depends upon the density ratio of 
ionizing photons to baryons as well as the shape of the 
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Fig. 9. — The fraction of carbon atoms in the C IV ionization 
state, as a function of redshift. Solid lines indicate the evolution 
for baryonic overdensities of p/p = 1 and p/p = 10, with the 
p/p = 10 curve peaking at z ~ 3.8 and the p/p = 1 curve peaking 
at z = 4.75. Dotted line indicates model for p/p = 10, but with 
additional local radiation from a galaxy at a distance of 100 kpc. 
The hatched region is excluded on the basis that it would lead to 
super-solar IGM abundances. 



background radiation spectrum. At high redshift, the 
ionizing background becomes less intense because of the 
declining QSO luminosity function, and softer because of 
an increasing contribution from galaxy light. Also, the 
average gas density increases as (1 + z) 3 , leading to a 
decline in the overall ionization level. 

The net effect of these changes is illustrated in Fig- 
ure which shows the evolution of fc rv from z ~ 
2 — 6. These calculations were performed using th e 
CLOUDY photoionization code ijFerland et al.1 [l998), 
with a background radiation spectrum containing both 



QSOs (/„ 



oc v 



) and galaxies (with 10% escape 



fraction of ionizing photons) after propegation through 
the IGM (|Haardt fc Madaul l200l|) . The spectra were 
normalized to be consistent with measureme nts of the 
proximity effect in QSOs ijScott et alJ [2000) and the 
distribution of flux transmission in the Ly-a forest 
ijMcDonald fc Miralda-Escudel20"o"l|) . The normalization 
becomes quite uncertain approaching z ~ 6; in this range 
I have estimated its value at by extrapolating from lower 
redshift data. The two solid lines indicate the C iv ion- 
ization fractions for gas with p/p = 1 and 10, where 
p = Qbp c (l + z) 3 . At z ~ 2, the intense radiation 
from QSOs ionizes most carbon to C V and other states 
higher than C iv. As one progresses through z ~ 4, 
the C IV fraction reaches its maximum and then declines 
rapidly at z > 4.5 as the carbon becomes more neutral. 

The total variation in fc iv spans 2-3 orders of magni- 
tude from z ~ 2 — 6 for a given overdensity. However, 
it approaches a maximum value of fc rv ~ 0.3 indepen- 
dent of density, though the corres ponding redshift does 
vary with p/p. As pointed out bv iSongailal l)2001t) . this 
fact can be used to set a lower limit on the overall carbon 
abundance at high redshift of Oc > flc iv/0.3 sa lx 10~ 7 . 
This can be compared to the overall baryon density 
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n b h 2 = 0.022 ijO'Meara et al.l200H) to obtain an volume- 
averaged carbon abundance: 




Here, mc and run represent the atomic mass of carbon 
and hydrogen, is the mean atomic mass in the IGM 
(assumed to be 1.3), and the last term represents the 
solar carbon abundance by num ber, assumed to be 3.3 x 
10" 4 lIGrevesse fc SauvaHfl998h . 

It should be heavily emphasized that this value does 
not represent a ubiquitous metallicity floor in the IGM. 
Rather it is a volume average over the entire redshift 
pathlength, which contains small regions with strong 
C IV absorption, and large regions with no detectable 
C IV . The expectation value of nc iv will be substan- 
tially lower than flc i V p c (l + z) 3 for most randomly cho- 
sen volumes, except for a few regions where it is much 
higher than this mean value. However, Equation does 
provide a global average which may be appropriate for 
comparing against the integrated metal content, for ex- 
ample, within an entire cosmological simulation volume. 

As a final note, it would also be interesting to measure 
the abundance of other carbon and/or silicon ions (e.g 
Si iv 1400, or C II 1334), since the ratios of these ions 
probe both relative heavy element abundances and the 
shape of the ionizing background spectrum (which is of 
particular interest during reionization at z > 6). Unfor- 
tunately the SNR of the present data are not sufficient 
to significantly constrain Si iv/C iv or other ion ratios, 
though this will be an important goal of future observa- 
tions. 

5.2. What does Qq iv Represent? 

The distribution of C IV shows several remarkable 
properties that might not be expected a priori. The most 
obvious of these is the near-constant value of fic iv from 
z = 2 to z = 6. The measurements are identical to 
within errors — even accounting for the 90% confidence 
intervals shown in Figure [H] f2c iv can vary by no more 
than a factor of ~ 3 over this complete redshift range. 
This comes despite the fact that fc rv can vary by up to 
3 orders of magnitude, and the expectation that cosmic 
abundances should decline on average towards higher z. 
In order to keep f2c iv constant, the evolution in ioniza- 
tion would need to precisely cancel an opposite evolution 
in the carbon abundance, a scenario that is possible, but 
which requires fine-tuning. 

Second, the C IV mass is distributed very unevenly. 
In the z ~ 5.72 - 6.21 bin, 43% of the C IV atoms are 
contained within a single system (at z = 5.823 — 5.830 
in SDSS1030). A similar segregation of C iv mass is 
seen at all redshifts. This is related to the fact that 
the CDDF slope is shallower than N^ v ; so that the 
integral flc iv °c / Ac rvf(Nc iv)dNc iv is dominated by 
the strongest system in any sample. 

At lower redshifts (z ~ 2.5 — 3.5), there is evidence 
that the strongest C iv absorption systems arise in the 
outer halos of star form ing galaxies ( jSimcoe et alJl2006i 
lAdelberger etld1l200l . They are often accompanied by 
heavy element absorption from many other ions including 
O vi , and their physical characteristics suggest a relation 



to energetic feedback from galaxy formation. Absorption 
from feedback-related syste ms comprises at least half of 
the f ig iv integral at z ~ 3 (iSimcoe et al.ll2006t iSongailal 
2006), and these lower redshift systems share some fea- 
tures with the z ~ 6 C iv systems. In particular the 
range of C iv column densities and complexity of velocity 
structure seem to be similar for the strongest C iv sys- 
tems at all epochs. 

This may provide one natural explanation for the con- 
stant value of f^c iv with redshift. If the C iv traced 
by the strong systems is produced locally by star form- 
ing galaxies, they can be enriched to 0.01 — 0.1-Z© lev- 
els on short timescales. In this case, ftc iv would be 
related more closely to the instantaneous metal produc- 
tion (i.e. star formation) rate at a given redshift, rather 
than the cumulative stellar mass formed prior to the 
epoch of observation. Much like flc iv, the star forma- 
tion rate dens ity appears to evolve fairly little between 
z = 2 - 6 (e.g.lGiavalisco et alJl200l iSteidel et al.lll999t 
iThomnson et al] 12006ft . Furthermore, if the C iv ab- 
sorbers reside near forming galaxies then the spectrum 
of ionizing radiation that they see would be augmented 
by local starlight, providing a buffer against large varia- 
tions in /civ 

In Figure [5J I illustrate this effect in the dotted curve, 
which shows the evolution of fc iv for the same Haardt 
& Madau diffuse background spectrum, but with an 
additional star forming galaxy located 100 kpc from 
the absorption system. The galaxy te mplate spectrum 
is tak en from the Starburst99 archive (Lei therer et alJ 
1999). It has a stellar age of 3 00 My r, and has been red- 
dened using a lCalzett i et alJ <l!994l) ex tinction law with 
E(B-V) = 0.155 <|Shaplev et al.l200l . The total stellar 
mass is 10 8 M Q . 

The galaxy has little effect on the ionization balance 
at z < 4.5, but at higher redshift the added galaxy flux 
does begin to contribute such that fc iv drops by only an 
order of magnitude from its peak, rather than 2—3 orders 
of magnitude as it would otherwise have done. It seems 
plausible that a decline in fc rv of a factor of < 10 could 
be coupled with a factor of < 10 decline in the carbon 
abundance at high redshift, resulting in weak observed 
evolution in Qc iv- 

The other possible scenario invokes a balance of sev- 
eral effects to produce a flat ttc iv(z), and does not re- 
quire a local origin for the observed C IV absorbers. In 
this interpretation, the C iv systems trace intergalac- 
tic regions of increasingly smaller overdensity at higher 
redshift, in essence picking out the density where the 
C iv ionization fraction (or fc rv times the heavy ele- 
ment abundance) peaks at any particular redshift. One 
can infer from a comparison of the solid lines in Figure 
[5] that this provides an independent mechanism to sup- 
press the variation in ionization parameter with redshift. 
It would require a fairly uniform distribution of carbon 
abundance across all densities of gas in the IGM, which 
although more difficult to achieve at high redshift, may 
still be possible. 

This explanation is difficult to test observationally, 
since at high redshift one also loses information about 
the local density field (because of saturation and blend- 
ing within the H I Lyman alpha forest). Naively one 
would expect d/V jdX for C iv to increase with redshift 
in this scenario (an effect which is not observed), un- 
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less the actual chemical abundances correlate with IGM 
density in just the right fashion. 

This "balancing" hypothesis convolves the evolu- 
tion of several competing factors (characteristic densi- 
ties, ionization parameters, and chemical enrichment), 
so its ultimate evaluation will be best accomplished 
through "observations" of numerical simulations which 
incorporate models of che mical feedback. Indeed, 
IQppenheimer fc Pavel l|2006|) have very recently run sim- 
ulations where galactic wind prescriptions motivated by 
local observations seem to strike the balance required to 
match the z — 2 — 5 observations. In their simulations 
there is evidence that flc w should start to turn down by 
z ~ 6 in contrast to the findings presented here. How- 
ever, the level of this discrepancy is probably not alarm- 
ing when one takes into account both our measurement 
uncertainties and uncertainties in the contribution from 
unresolved small-scale C iv structures in the simulations. 

6. CONCLUSIONS 

This paper presents moderate resolution spectra of two 
z ~ 6 QSOs, taken with GNIRS on Gemini South. These 
spectra were systematically searched for high redshift in- 
tergalactic C iv absorption, using objective selection cri- 
teria. Comprehensive testing was performed to quantify 
the C IV sample incompleteness as a function of Nq iv , as 
well as the rate of false positive C iv identifications. Over 
the primary redshift range of z = 5.72 — 6.21, the sam- 
ple is largely complete at Nq i V ~ 10 14 cm -2 , and drops 
uniformly to about 5 — 10% at Nq iv ~ 10 13 ; at lower 
redshift the completeness is somewhat worse because the 
data have a lower SNR. However, this sensitivity is ade- 
quate to construct a rough estimate of the C iv column 
density distribution f(Ncrv), and also to measure the 
integrated density Qc iv at a level two to three times be- 
low the observed value at z = 2 — 5. The main results 
are as follows: 

1. Nine C iv absorption lines are identified between 
z = 5.72 and z = 6.21. Three of the these lines 
are contained within a single strong complex at 
z = 5.827 ± 0.004 toward SDSS1030+0524. This 
single system contains 43% of the C iv atoms along 
the total combined redshift path. Three additional 
C iv lines are detected in a lower redshift bin rang- 
ing from z — 5.40 — 5.72. After adjusting for ab- 
sorption pathlength and incompleteness, the den- 
sity of lines is similar in both redshift bins. 

2. The spectral region at z > 6 in SDSS1030+0524 
is statistically consistent with a complete lack of 
C iv absorption. However, with only one sightline 
sampled at present this could easily be the result 
of a statistical fluctuation. 

3. A coarse-binned representation of the z ~ 6 
C iv column density distribution is consistent with 
the redshift-in variant form rep orted over the z = 
2-5 range bv lSongailal l|2001lK 

4. An extension of the ficiv(^) plot to z ~ 6 does 
not show evidence of a downturn in the integrated 
C iv abundance, at least within the factor of ~ 3 
allowed by measurement errors. 



Over such a large redshift range, it is remarkable to 
derive a constant form for the C iv column density distri- 
bution and also for Qc iv The baryonic density, ionizing 
background spectrum, and metallicity should all evolve 
substantially between z = 2 — 6, so it would represent 
an interesting coincidence if these effects all cancelled. 
In principle it is possible to mask the effect of evolu- 
tion in the C IV ionization fraction if the C iv absorbers 
are associated with regions of lower overdensity at higher 
redshifts. In this case one might anticipate an increase 
in dAf jdX at higher z, but this is not observed. Because 
of the complex interplay of many evolving factors, a full 
evaluation of this explanation for the constancy of fic iv 
with z requires comparison with numerical simulations 
(e.g. IQppenheimer fc Davel l2006). 

An attractive alternative explanation is that the Qq w 
measurements are dominated by strong absorbers whose 
properties are affected by local chemical and radiative 
feedback from nearby galaxies. This association appears 
to be present at z ~ 2.5, where luminous galaxies are 
seen directly in the neighborhood of many/most strong 
C iv systems. Although the same galaxies cannot be 
seen at high redshift, they would in principle provide a 
mechanism to rapidly enrich the C iv absorbing regions, 
and their radiation would also regulate the ionization 
balance of the carbon, preventing it from becoming too 
neutral for detection in C IV . 

Local feedback models become more compelling at the 
highest redshifts, where the elapsed time betewen the 
onset of reionization (z ~ 11) and the epoch of observa- 
tion becomes quite short. Over timescales of 400 — 500 
Myr, there would be just enough time to form an early 
galaxy population, build up its stellar mass, and eject 
heavy elements from its earliest supernovae into a region 
< 100 kpc in size. In this scenario there could be large 
metal-free regions of the IGM, but C IV observations at 
present sensitivities would only probe the regions where 
local production was important. In time, IR spectra with 
higher resolution and SNR should be able to track the 
evolution of weaker C iv lines in the tenuous IGM, and 
these systems will provide a more sensitive test of local 
versus gloal enrichment scenarios in the early universe. 
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TABLE 1 

C iv Lines Detected Toward SDSS1306+0356 



ZC iv 
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aw 
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EW Ratio 
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log(AT c IV ) 
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TABLE 2 

C iv Lines Detected Toward SDSS1030+0524 



ZC iv 


a z 
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TABLE 3 

C iv Column Density Distribution by Redshift 



Az 




f 2 

J comp 


13.00 

AX cff 3 


- 13.75 
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log(AT c 
log/(iVciv) 4 | 
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/comp 
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- 14.50 
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log/(JV C iv) 
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2.98 
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0.412 
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-13.27 


0.624,0.345 


1.440 
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-14.27 


5.72-6.21 


3.35 


0.500,0.318 


1.247 
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-13.16 


0.880,0.709 


2.546 
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-14.34 


5.90-6.21 


1.64 


0.379,0.350 


0.578 


4 


-12.83 
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< -14.48 



1 Full absorption path (denned by Equation^ along redshift range assuming 100% sample completeness. 
2 Compleness along each sightline for given range in redshift and Nc iv- First value is for SDSS1306+0356, 
second is for SDSS1030+0524. 
3 Effective redshift path for bin, determined by applying the corrections given in previous column. 
4 Column density distribution function value. 



